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Abstract 

For special coupling ratios, the one-dimensional (ID) s = 1/2 Heisenberg 
model with antiferromagnetic nearest and next-nearest neighbor interactions 
has a pure dimer ground state, and the ID s = 1 Heisenberg model with 
antiferromagnetic bilinear and biquadratic interactions has an exact valence- 
bond-solid ground state. The recursion method is used to calculate the T = 
spin dynamic structure factor for both models and, for the s = 1/2 model, 
also the dimer dynamic structure factor. New results for line shapes and 
dynamically relevant dispersions are obtained. 
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Correlated quantum fluctuations in the ground state are a generic feature of quantum 
many-body systems. They make it hard to take finite-size effects into account in compu- 
tational studies of zero-temperature dynamical properties. Interestingly, there are several 
known cases, where a relatively simple-structured ground state is stabilized by competing 
terms in the microscopic Hamiltonian. This ground state may or may not be long-range 
ordered. The essential attribute is that its fluctuations are not correlated or only over a 
short distance on the lattice. This phenomenon typically occurs with no accompanying sim- 
plification in the excitation spectrum or any dynamical quantity. Nevertheless, any such 
situation provides an unsuspected window for dynamical studies which promise to be much 
less plagued with finite-size effects than is typically the case. The recursion methodB in 
combination with recently developed techniques of continued-fraction analysis^ is an ideal 
calculational tool for that purpose, the key property being that it extracts the dynamical 
information from the ground-state wave function. 

In a previous paper we have reported the study of one such case, namely the one- 
dimensional (ID) spin-s XYZ model.il In a magnetic field of particular strength, this model 
has a product ground state with spontaneous ferro- or antiferromagnetic long-range order 
perpendicular to the field. Here we present new results for two different models with simple- 
structured ground states. 

The first model is the ID s = 1/2 Heisenberg antiferromagnet with competing nearest 
and next-nearest neighbor interactions, 

N 

H = J2ViSi ■ Siu + J2S1 ■ S l+2 }, (1) 
1=1 

with an even number of spins and periodic boundary conditions. This system undergoes a 
T = phase transition at J^j J\ — 0.25 from a spin-fluid phase to a phase with spontaneous 
dimer long-range order.i In the (critical) spin-fluid phase, the correlations of the quantum 
fluctuations are particularly strong. In the dimer phase, their continued presence, albeit 
much attenuated, manifests itself, for example, in the finite-size splitting of the ground- 
state doublet. The exception is the special coupling ratio J2/J1 = 0.5, where the pure 



dimer ground-state is realized. Bo Here the ground-state energy per site is size-independent: 
E /N=-lJ 1 . 

The two (translationally invariant) dimer ground-state wave functions can be expressed 
in terms of products of singlets formed by pairs of nearest neighbor spins: 

l*±> = (2 ± (-l/2)^)~ V2 {!<&!> ± |$ 2 >}, (2) 

where 1$!) = [1, 2] [3, 4] • • • [N-l, N], |$ 2 ) = [2, 3] [4, 5] • • • [N, 1], [1, 1+1} = {| U)-| tf)}/y/2. 

The dimer order parameter, D = iV _1 £;(— l) l Di, D\ = SfS^ + SfS^, has a nonzero 
expectation value, (D) = ±1/2 in the (non-orthogonal) symmetry-breaking states |$i) and 
| $2)- The order-parameter correlation function in this case is not a two-spin correlation 
function, {(SfSf+n) = for \n\ > 1), but a four-spin correlation function: (DiD i+n ) = 
(— l) n /4 for n ytz 0. Hence it will be instructive to compare the spin dynamic structure 

factor S zz (q,u) and the dimer dynamic structure factor Sdd(<?, u>), i.e. the function 

+00 

S AA (q^)= J dte^(A q (t)A_ q ) , (3) 

—00 

where A q stands for the spin fluctuation operator, S q = N~ l l 2 Y>ie %ql Sf or the dimer fluctu- 
ation operator, D q = N^ 2 H l e iql [D l - (D,)]. 

By means of the recursion method^ in combination with a strong-coupling continued- 
fraction analysis,^ we calculate the dynamically relevant excitation spectra and the spectral- 
weight distributions of these two functions. The recursion algorithm in the present context 
is based on an orthogonal expansion of the wave function \^ q (t)) = A q (—t)\$). It produces 
(after some intermediate steps) a sequence of continued-fraction coefficients Af(q), Af(q), . . . 
for the relaxation function, 



Z+ ... 

which is the Laplace transform of the symmetrized correlation function §l(A q (t)A- q ) / 
(A q A_ q ). The T = dynamic structure factor (|3|) is then obtained from (|j) via 



S AA (q,u) = 4(A^- q )e(u)^^c* A (q,e-iw)] . (5) 

The simple iV-dependence of the dimer ground-state wave functions (|2|) offers the advantage 
that we can compute a significant number of iV-independent coefficients A^(g). These data 
are the input to the well-tested strong-coupling continued-fraction reconstruction.!'! 

The results for the frequency-dependence of the dynamic structure factor S zz (q,u>) at 
the wave numbers realized for N = 18 are displayed in Fig. |l[ The set of curves is perfectly 
compatible with a function S zz (q,uj) that varies smoothly in q as well as in u. For every 
q- value we observe a single peak in the frequency range of interest. This peak is very broad 
at q near zero or n. The width shrinks as q approaches 7r/2 from either side. 

At q = 7r/2 (not realized for N = 18) the dynamically relevant excitation spectrum 
reduces to a single mode. The states S*, 2 \$>±) are, in fact, known to be exact eigenstates of 
the system.!'! The dynamically relevant dispersion of S zz (q,u) is symmetric about q = ir/2, 
where it has a smooth maximum at frequency u/J\ = 1.0. For q = and q = 7r it has 
smooth minima at u/J\ c± 0.5. 

On the basis of a variational calculation for the pure dimer state, Shastry and Sutherland© 
obtained an excitation spectrum for this model which consists of a continuum of two-defect 
scattering states with a lower boundary e(q) — J\{\ — | cosg|) and, for the restricted range 
0.367T < q < 0.647T of wave numbers, a branch of defect bound states, which emerges from 
the lower continuum boundary and has a smooth maximum reaching up to u)/J\ = 1 at 
q = tt/2. 

In this context, our results suggest that for q near zero or ir, the spectral weight of 
S zz (q, ^) is distributed over a broad frequency range of two-defect scattering states. As q 
approaches 7r/2 from either side, the spectral weight is shifted gradually to the two-defect 
bound state. 

Our result for the dynamic structure factor SDd(q,uj) is plotted in Fig. |] for the same 
set of wave numbers. The line shapes and peak positions resemble those shown in Fig. [I] 
for S zz (q, lj), but there are some notable differences: The spectral weight in Snn^q.uj) is 
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concentrated at somewhat higher energies. The intensity at q = is nonzero. The shape of 
the dimer dispersion is different. S , dd( 7I "/2, uj) does not reduce to a single line. In Sdd(k,u) 
the spectral weight is shared between a continuum and a 5-peak at u — 0.0 The latter 
contribution reflects the presence of dimer long-range order in the ground state. 

The Hamiltonian of the ID s = 1 model with isotropic bilinear and biquadratic exchange, 
our second example, is most conveniently expressed in the form 

N 

H-y = J ^{cos7S/ • S; + i + sin7(Si • S m ) 2 } (6) 
i=i 

with a single parameter —it < 7 < n. More than a decade of research on this model has 
established a T = phase diagram consisting of the short-range ordered Haldane phase 
and three phases with dimer, trimer, and ferromagnetic long-range order.0 In the Haldane 
phase, which includes the Heisenberg antiferromagnet (7 = 0), the ground state is a non- 
degenerate singlet state (St = 0) separated by a gap from the threshold of the excitation 
spectrum. 

At the parameter value 7 = arctan(l/3) ~ 18.4° within this phase, the ground-state wave 
function is exactly known Jil It is a realization of the so-called valence-bond-solid (VBS) wave 
function, which can be assembled from the same parts as the dimer state (Q). The spin 1 at 
each lattice site is expressed as a spin- 1/2 pair in a triplet state. The singlet-pair forming 
valence bond involves one fictitious spin 1/2 from each of two neighboring lattice sites. The 
VBS state can then be regarded as a chain of valence bonds linking successive symmetrized 
spin-1/2 pairs, which are given by the S T = vector of the triplet on each lattice site. 

The static spin correlation function in the VBS state,@ (SfSf +n ) = |(-l) n 3 _ l n l, (n ^ 0), 
reflects magnetic short-range order with a very short correlation length: £ = l/ln3 ~ 0.91. 
The static structure factor is then a non-singular function with a smooth minimum at q = 
and a smooth maximum at q = tt: 

S zz (q) =2(l-cosg)/(5 + 3cosg). (7) 

The simple structure of the VBS state makes this quantity free of finite-size effects for 
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q = 2irl/N, I = 0, . . . , iV — 1. Again, this simplification does not extend to the excitation 
spectrum and the dynamical properties. 

In Fig. H we have plotted S zz (q, to) as obtained from a strong-coupling continued- fraction 
analysis with coefficients extracted from the N = 12 VBS wave function. At each value of q 
the spectral weight of S zz (q, uo) is found to be concentrated in a single peak with symmetric 
line shape. The peak frequency decreases monotonically with q. The suggested gap value 
at q = tt is AE/J ~ 0.66. The linewidth tends to be very small at q near tt, where the peak 
frequency is lowest and the intensity highest. It gains considerably in breadth at q near 0, 
where the peak frequency is higher and the intensity much lower. 

The monotonic g-dependence of the dynamically relevant dispersion for S zz (q,u) in the 
VBS state (7 ~ 18.4°) is markedly different from the corresponding quantity in the Heisen- 
berg antiferromagnet (7 = 0), which belongs to the same phase. In the Heisenberg case, 
the dispersion has a smooth maximum at q ~ tt /2 and smooth minima of unequal height at 
q = and q = 7T.E1 

This work was supported by NSF Grant DMR-93-12252 and by the NCSA at Urbana- 
Champaign. 
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FIGURES 

FIG. 1. Dynamic structure factor S zz (q,uj) vs uj for q = 2nl/N, I = 0,1,..., N/2 in the 
dimer ground state |$+) of the Hamiltonian (|l|) at J2 = 0.5 and J\ = 1 with N = 18, obtained 
via strong-coupling continued-fraction reconstruction based on the coefficients Ai,...,Ag and a 
Gaussian terminator as explained in Refs. |],|7]. 

FIG. 2. Dynamic structure factor Sdd(q,w) vs uj for q = 2nl/N, I = 0, 1, . . . , N/2 in the 
dimer ground state |3>+) of the Hamiltonian (|l|) at J2 = 0.5 and J\ = 1 with iV = 18, obtained 
via strong-coupling continued-fraction reconstruction based on the coefficients Ai,...,Ag and a 
Gaussian terminator. 

FIG. 3. Dynamic structure factor S zz (q,u) vs u> for q = 2irl/N, I = 0, 1, . . . , iV/2 in the VBS 
ground state of the model system (P) at J = 1 and 7 = arctan(l/3) with N = 12, obtained via 
a strong-coupling continued- fraction analysis based on the coefficients A\,...,Aq and a Gaussian 
terminator 
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